The ISME Journal (2014) 8, 1153-1165 

© 2014 International Society for Microbial Ecology All rights reserved 1751-7362/14 




www.nature.com/ismej 

ORIGINAL ARTICLE 

Endospores of thermophilic bacteria as tracers 
of microbial dispersal by ocean currents 

Albert Leopold Miiller 1 ' 2 , Julia Rosa de Rezende 3 ' 4 , Casey RJ Hubert 4 , Kasper Urup 
Kjeldsen 3 , Ilias Lagkouvardos 1 , David Berry 1 , Bo Barker j0rgensen 3 and Alexander Loy 1 ' 2 
^Division of Microbial Ecology, Department of Microbiology and Ecosystem Science, Faculty of Life Sciences, 
University of Vienna, Vienna, Austria; 2 Austrian Polar Research Institute, Vienna, Austria; ^Center for 
Geomicrobiology Department of Bioscience, Aarhus University, Aarhus, Denmark and 4 School of Civil 
Engineering and Geosciences, Newcastle University, Newcastle Upon Tyne, UK 



Microbial biogeography is influenced by the combined effects of passive dispersal and 
environmental selection, but the contribution of either factor can be difficult to discern. As 
thermophilic bacteria cannot grow in the cold seabed, their inactive spores are not subject to 
environmental selection. We therefore conducted a global experimental survey using thermophilic 
endospores that are passively deposited by sedimentation to the cold seafloor as tracers to study 
the effect of dispersal by ocean currents on the biogeography of marine microorganisms. Our 
analysis of 81 different marine sediments from around the world identified 146 species-level 16S 
rRNA phylotypes of endospore-forming, thermophilic Firmicutes. Phylotypes showed various 
patterns of spatial distribution in the world oceans and were dispersal-limited to different degrees. 
Co-occurrence of several phylotypes in locations separated by great distances (west of Svalbard, 
the Baltic Sea and the Gulf of California) demonstrated a widespread but not ubiquitous distribution. 
In contrast, Arctic regions with water masses that are relatively isolated from global ocean 
circulation (Baffin Bay and east of Svalbard) were characterized by low phylotype richness and 
different compositions of phylotypes. The observed distribution pattern of thermophilic endospores 
in marine sediments suggests that the impact of passive dispersal on marine microbial 
biogeography is controlled by the connectivity of local water masses to ocean circulation. 
The ISME Journal (2014) 8, 1153-1165; doi:10.1038/ismej.201 3.225; published online 19 December 2013 
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Introduction 

Microorganisms display distinct biogeographic pat- 
terns (reviewed in Foissner, 2006; Green and 
Bohannan, 2006; Martiny et al., 2006; Ramette and 
Tiedje, 2007; Lindstrom and Langenheder, 2012), yet 
the mechanisms controlling their distribution in the 
environment are difficult to distinguish and thus not 
well understood. Four fundamental processes — 
selection, drift, dispersal and mutation — have been 
proposed for creating and maintaining microbial 
biogeographic patterns (Hanson et al, 2012), updat- 
ing the classical concept of dispersal, speciation and 
extinction as the main factors determining biogeo- 
graphy Regardless of theoretical framework, 
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dispersal of microbial cells has a central role in shaping 
the spatial distribution of microbial biodiversity 
(Green and Bohannan, 2006; Fierer, 2008). In its 
strictest sense, microbial dispersal is defined as the 
physical movement of cells between two locations, 
but an extended definition additionally includes 
successful establishment — that is, physiological 
activity and growth of migrated cells — at the receiv- 
ing location (Hanson et al., 2012). The existence of 
physical dispersal barriers for microorganisms has 
traditionally been questioned. Under the Baas 
Becking paradigm of microbial cosmopolitanism — 
'everything is everywhere, but, the environment 
selects' — it is hypothesized that microorganisms 
possess unlimited dispersal capabilities due to their 
large population sizes and short generation times. 
In this paradigm, environmental factors are the sole 
determinants of observed microbial distribution 
patterns (Baas Becking, 1934; Finlay, 2002; Fenchel 
and Finlay, 2004). More recently, multiple studies 
have put forward evidence for dispersal limitation 
among microorganisms (Papke et al., 2003; Whitaker 
et al., 2003; Green et al., 2004; Reche et al., 2005; 
Martiny et al., 2006; Ghiglione et al., 2012; Sul et al., 
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2013). Controversy remains, however, partly due to 
differences in definitions of dispersal and the 
operational taxonomic units used to classify micro- 
organisms (Hanson et al., 2012), but mostly because 
the biogeography of free-living microorganisms is 
the result of a combination of various evolutionary 
and ecological processes that make it difficult to 
unravel the relative influence of passive transport 
(Martiny et al, 2006; Ramette and Tiedje, 2007). We 
use the term dispersal here to describe physical 
movement by passive transport, but not colonization 
of the new location. 

The presence of dormant endospores of thermo- 
philic members of the bacterial phylum Firmicutes 
in cold marine sediments is the result of passive 
transport from warm environments and has been 
proposed as a natural model for selectively investi- 
gating the dispersal of microbial cells in the oceans 
(Bartholomew and Paik, 1966; Isaksen et al., 1994; 
Hubert et al., 2009; de Rezende et al., 2013). The 
inactivity of these model organisms enables their 
biogeography to be investigated largely without any 
confounding influence of environmental selection. 
Therefore, any observed non-random spatial distri- 
bution patterns of these organisms should be 
directly attributable to the influence of dispersal 
limitation. Ocean currents and eventual sedimenta- 
tion have been invoked as dispersal vectors for 
supplying thermophilic Firmicutes spores to Arctic 
fjord sediments off the coast of Svalbard at an 
estimated rate exceeding 10 8 spores per square meter 
per year (Hubert et al., 2009). Once deposited in the 
cold sediment, these spores lie dormant as members 
of the rare biosphere and the 'microbial seed bank' 
(Pedros-Alio, 2012; Gibbons et al, 2013), but they 
can be induced to germinate rapidly during sedi- 
ment incubation experiments at high temperature 
(Hubert et al, 2009). 

Enrichment of thermophilic anaerobes from dor- 
mant spores in high-temperature incubation experi- 
ments enables selective focus on a specific, 
culturable part of the rare biosphere that is of 
exceptional interest for studying long-term and 
long-distance dispersal. Bacterial endospore disper- 
sal represents an upper boundary for the dispersal 
capabilities of vegetative cells, as bacterial endo- 
spores are metabolically inert, highly stress resistant 
and able to survive unfavorable conditions for long 
periods (Nicholson et al., 2000). For example, a half- 
life of up to 440 years was estimated for endospores 
of thermophilic sulfate-reducing bacteria deposited 
in Aarhus Bay in the Baltic Sea (de Rezende et al., 
2013); these viable endospores decreased in abun- 
dance with depth but were still detected at 6.5 m 
below seafloor in 4500-year-old sediment, suggest- 
ing a life span that could allow the global dispersal 
of spores through thermohaline circulation that 
fully connects the world oceans on a time scale of 
-1000-2700 years (DeVries and Primeau, 2011). 
Diverse populations of dormant bacterial spores 
become active when sediment is heated (to 50-60 °C). 



Once germinated, the vegetative cells catalyze the 
mineralization of organic matter via extracellular 
enzymatic hydrolysis, fermentation and sulfate 
reduction (Hubert et al., 2010). So far, the phyloge- 
netic composition of these dormant, thermophilic 
communities has only been studied in a few 
locations, that is, the sediments of West Svalbard 
fjords and Aarhus Bay. Even though these locations 
are ~ 3000 km apart, they share at least two thermo- 
philic Desulfotomaculum phylotypes with identical 
16S rRNA and dissimilatory (bi)sulfite reductase 
[dsrAB) gene sequences (de Rezende et al., 2013). 
Although this is the first intriguing evidence for 
long-distance dispersal of these thermophilic endo- 
spores, their large-scale biogeography in the global 
ocean remains unexplored. 

We thus analyzed the richness, phylogeny and 
distribution of endospores of thermophilic bacteria 
in marine sediments on global and regional scales to 
address the following questions: are thermophilic 
spores randomly distributed, which would be 
indicative of unlimited dispersal, or does the 
biogeography of thermophilic spores show signs 
of dispersal limitation? If these spores are not 
randomly distributed, how do local hydrography, 
sedimentation and major ocean currents impact 
their distribution? To investigate the biogeography 
of thermophilic spores, sediment samples from 81 
locations around the world ocean, including regio- 
nal Arctic samples from the Svalbard archipelago 
and the Baffin Bay that connect differently to global 
ocean circulation, were amended with organic 
substrates, pasteurized and incubated at 50 °C under 
anoxic conditions. Pasteurization kills the vegeta- 
tive community of mesophilic or psychrophilic 
microorganisms. Shifts in microbial community 
activity and composition, owing to germination 
and growth of thermophilic endospores during the 
incubation, were then monitored by measuring 
sulfate reduction rates and by using multiplex 
pyrosequencing of bacterial 16S rRNA gene amplicons. 
Subsequently, the phylogeny and biogeography of 
enriched thermophilic phylotypes was analyzed. 



Materials and methods 

Marine sediment samples 

The sample set comprised marine sediments from 
81 locations around the world ocean, including two 
regional Arctic sample sets from Svalbard fjords and 
the Baffin Bay and samples from hydrothermally 
influenced sediments in the Guaymas Basin 
(43-150 °C, Gulf of California) that comprises 
potential source environments of endospore-forming 
thermophiles (Figure 1, Supplementary Table Si). 
Most sediment samples were collected from the 
seafloor surface (0-10 cm below seafloor) at coastal 
or deep sea, open ocean sites with in situ temperatures 
ranging from 0 to 30 °C. Samples were stored at 4 °C 
or frozen at — 20 °C until germination experiments. 
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Figure 1 Global and regional maps show the sediment-sampling sites and selected major ocean currents. Circle sizes represent richness 
of thermophilic Firmicutes endospore phylotypes (crosses indicate a richness of 0). Symbol color indicates whether thermophilic sulfate 
reduction was detected during the high temperature incubation (green = positive, red = negative). Global map (top) shows the global 
thermohaline circulation (warm/surface currents in red, cold/deep water currents in blue; adapted and simplified from Rahmstorf 
(2002)). Broad geographic regions studied are highlighted in lighter blue (WS, West Svalbard; ES, East Svalbard; BB, Baffin Bay; NE, 
Northern Europe; GB, Guaymas Basin; SA, South Asia; SP, South Pacific). Map of Baffin Bay (bottom left): local currents reproduced from 
Lloyd et al. (2005); BC, Baffin Island Current; WGC, West Greenland Current; LC, Labrador Current. Map of Svalbard (bottom right): local 
currents reproduced from Knies et al. (2007); WSC, West Spitsbergen Current; ESC, East Spitsbergen Current. 
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Sediment incubation and sulfate reduction 
measurements 

Sediment was homogenized and mixed with sterile 
anoxic synthetic seawater medium at a 1:2 (w/w) 
ratio under constant flow of N 2 gas. The sediment 
slurries were amended with organic substrates: 
formate, lactate, acetate, succinate, propionate, 
butyrate, ethanol (each to a final concentration of 
0.5 mM) and/or with freeze-dried Spirulina cells 
(1.5 gl" 1 ). Two 12-ml aliquots of slurried sediment 
were transferred to Hungate tubes under constant 
flow of N 2 and pasteurized for 20min at 80 °C and 
incubated in parallel at 50 °C. Before the incubation, 
one aliquot received ~720kBq 35 S-labeled carrier- 
free sulfate tracer for sulfate reduction measure- 
ment. The incubations were sub-sampled by 
syringe-needle after 0, 56, 72 and/or 120 h. Aliquots 
of 3 ml from the 35 S-sulfate slurries were mixed with 
6 ml 20% zinc acetate solution and stored at - 20 °C 
until sulfate reduction was determined using a 
single-step cold chromium distillation method 
(Kallmeyer et ah, 2004). From the non-radioactive 
incubations, 2 ml subsamples were pelleted by 
centrifugation and frozen at - 20 °C for subsequent 
DNA extraction. 



Multiplex amplicon pyrosequencing 
DNA was extracted from sediment slurries using the 
PowerSoil DNA isolation kit (MO BIO Laboratories, 
Inc., Carlsbad, CA, USA). Bacterial 16S rRNA gene 
amplicons were obtained and supplied with bar- 
codes and pyrosequencing adaptors using a two-step 
PCR approach with low cycle numbers (20 + 5) and 
pooling of triplicate PCR reactions to reduce 
variability associated with barcoded pyrosequen- 
cing primers (Berry et ah, 2011). The PCR primers 
targeting most bacteria (909F: 5'-ACTCAAAKGAAT 
WGACGG-3 / , 1492R: 5 / -NTACCTTGTTACGACT-3 / ) 
and conditions were described previously (Berry 
et ah, 2011). Final PCR products were purified using 
the Agencourt Ampure XP system (Beckman Coul- 
ter, Vienna, Austria) and the DNA concentration 
was determined using a Quant-iT PicoGreen 
dsDNA Assay (Invitrogen, Darmstadt, Germany). 
For better coverage of samples with higher 
expected bacterial diversity, amplicons were then 
pooled at a 2:1 ratio of Oh and 120 h time points, 
respectively. Sequencing was performed on a GS 
FLX or GS FLX+ instrument using Titanium 
chemistry (Roche, Mannheim, Germany) by the 
Norwegian High-Throughput Sequencing Centre 
(Oslo, Norway) or by Eurofins MWG Operon 
(Ebersberg, Germany). In total, 1 493 577 reads with 
an average length of 567 nucleotides (nt) were 
received. Sequences were trimmed and erroneous 
sequencing reads removed using the PyroNoise 
implementation in mothur (Quince et ah, 2009; 
Schloss et al., 2009) and sorted according to 
barcode using QIIME (Caporaso et al., 2010), 
yielding a total of 1196 847 usable reads with an 



average length of 353 nt. A 97% identity threshold 
was used for clustering reads into phylotypes with 
UCLUST (Edgar, 2010). Representative sequences 
were aligned with mothur using the Needleman- 
Wunsch pairwise alignment method default 
settings (Schloss et al., 2009). Chimeras were 
detected using Chimera Slayer (Haas et al., 2011) 
and excluded from further analysis. Pyrosequen- 
cing data are archived at the NCBI Sequence Read 
Archive under accession SRP028774. 



Identification of putative thermophilic endospore 
phylotypes 

Two criteria were used to identify putative thermo- 
philic endospore phylotypes. First, species-level 
phylotypes had to be significantly enriched in at 
least one sediment sample after incubation at high 
temperature. This criterion is based on the assump- 
tion that thermophilic endospores will survive the 
initial pasteurization and germinate and grow 
during the ensuing 50 °C incubation. Significant 
enrichment of phylotypes was determined using a 
two-proportion T-test and P-values were corrected 
for multiple comparisons using the false discovery 
rate method in R (Benjamini and Hochberg, 1995) to 
account for uncertainty due to sequence sampling 
depth. Corrected P-values < 0.01 were considered 
significant. Second, phylotypes had to be affiliated 
with the phylum Firmicutes, to which all known 
endospore-forming bacteria belong (it was shown 
recently that endosporulation likely evolved only 
once at the base of the Firmicutes tree) (Abecasis 
et al., 2013). This approach does not survey all 
inactive cells in the sampled environments, but 
hones in on a physiological subset of the rare 
biosphere (Lennon and Jones, 2011). The clear 
advantage of this strategy is that taxa that were 
active and subject to environmental selection in situ 
are excluded from the analysis, allowing passive 
dispersal to be evaluated in isolation. Representa- 
tive sequences of enriched Firmicutes phylotypes 
were then automatically aligned using the web- 
based SINA aligner (Pruesse et al, 2012) and 
imported into the ARB-SILVA database SSU Ref 
NR 111 (Quast et al., 2013) using the ARB software 
package (Ludwig et al., 2004). The alignment was 
manually curated and used to re-cluster the repre- 
sentative sequences into species-level phylotypes of 
^97% sequence similarity using the average neigh- 
bor algorithm in mothur. A phylotype was called 
present at a location if at least one sequence of this 
phylotype was detected before and/or after the 120 h 
incubation of sediment from this location. 



Phylogeny 

A maximum likelihood (RAxML) tree was calcu- 
lated with almost full-length 16S rRNA sequences 
(^1400 nt) of closely related reference bacteria or 
environmental clones based on 1222 alignment 
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positions by using a 50% sequence conservation 
filter for bacteria. Using the ARB Parsimony 
Interactive tool, the short amplicon pyrose- 
quences were subsequently added to this tree 
one at a time by using the 50% sequence 
conservation filter and alignment filters covering 
the individual length of each representative 
phylotype sequence, without changing the overall 
tree topology. Trees were visualized using iTOL 
(Letunic and Bork, 2006). 



Network analysis 

Two different kinds of networks were built based 
on the phylotype presence/absence matrix. First, a 
network of phylotype co-occurrence was produced 
for phylotypes present in at least five sites and 
with a minimum Spearman correlation coefficient 
of 0.6 (and P< 0.0001 based on permutation 
testing) (Barberan et al., 2012). Focus on phylo- 
types that occurred at multiple sites reduced 
network complexity and facilitated identification 
of the core thermophilic endospore community. 
Second, a network of sites based on phylotype 
diversity was produced based on Bray-Curtis 
similarity of at least 0.6. Networks were plotted 
using the 'network' package in R (Butts et al., 
2012). 



Distance-decay analysis 

Jaccard similarities were calculated from the 
phylotype presence/absence matrix using the 
Vegan' package in R (Oksanen et al., 2012). 
Distances were calculated from latitude and long- 
itude coordinates of sampling sites using the 
Vincenty inverse formula for ellipsoids and World 
Geodetic System (WGS) 84 ellipsoid parameters 
and calculated using R (R Development Core 
Team, 2008). A linear model was fit to the data 
and the statistical significance of the resulting 
regression parameters was evaluated by analysis of 
variance in R. To determine which phylotypes 
were dispersal-limited, a permutation-testing 
approach was used. The geographical dispersal 
of each phylotype was calculated as the mean 
distance between sites at which it was found. 
A null distribution of mean distances was produced 
by randomly selecting from all sites an equally 
sized subset (e.g. if a phylotype was found at seven 
sites, then seven sites would be randomly 
selected). The mean distance between the random 
subsets was calculated, and this process was 
repeated for 10 000 re-samplings. The probability 
that the observed mean distance was due to chance 
was calculated and corrected for multiple compa- 
risons and P-values <0.05 were considered 
significant. 

Geographical maps were drawn using Generic 
Mapping Tools (Wessel and Smith, 1998). 



Results and discussion 

Non-uniform distribution of thermophilic Firmicutes 
endospores in marine sediments 

We studied the biogeography of thermophilic 
endospores in Arctic and other permanently cold 
or temperate marine sediments, because longevity 
of endospores makes them ideal study objects 
for understanding the time-averaged impact of 
dispersal on marine microbial biogeography. 
While there might be differences in endospore 
survival between different bacteria (Nicholson 
et al., 2000; McKenney et al., 2013), results are 
interpreted under the assumption that differences 
in endospore distribution are dependent on 
dispersal and not on environmental selection. 

Anoxic, high-temperature incubations of pasteur- 
ized marine sediment (mostly from the top 
few centimeters) (Supplementary Results and 
Discussion, Supplementary Figure Si) from 81 
different locations (Supplementary Table Si) led 
generally to a reduction in bacterial alpha diversity 
due to germination and growth of thermophilic 
endospores combined with the death and DNA 
decay of vegetative cells (Supplementary Materials 
and Methods, Supplementary Table S2). Principle 
coordinate analysis of weighted UniFrac distances 
confirmed a shift in phylogenetic composition 
of bacterial 16S rRNA genes after the incubation of 
most, but not all, sediment samples (Supplementary 
Figure S2). In total, we identified 146 thermophilic 
endospore phylotypes (hereafter called 'thermo- 
spore phylotypes') across all high-temperature 
incubations (Supplementary Table S3). Thermo- 
spore phylotypes were detected in samples from 
almost all of the investigated locations [n = 79/81) 
(Figure 1). Most of the 146 thermospore phylotypes 
were affiliated with the orders Clostridiales 
(61.0%) and Bacillales (17.8%); the most repre- 
sented families were Clostridiaceae (32.2%), 
Peptococcaceae (17.8%) and Bacillaceae (14.4%) 
(Supplementary Figure S3, Supplementary Table S3). 
The phylogenetic identity and potential physiology 
of thermospore phylotypes corroborate and greatly 
expand upon previous studies that investigated 
thermophilic endospores in marine sediments 
(Bartholomew and Paik, 1966; Isaksen et al, 1994; 
Vandieken et al, 2006; Hubert et al, 2009, 2010; Ji 
et al., 2012; de Rezende et al., 2013) (Supplementary 
Results and Discussion). A screening of environmental 
amplicon databases showed that 16S rRNA gene 
sequences of the identified thermospore phylotypes 
are very rarely detected by environmental surveys in 
the marine environment (Supplementary Results and 
Discussion, Supplementary Table S4), which is likely 
consistent with the widespread occurrence of 
endospore taxa in low relative abundance in many 
habitats. 

The different sampling locations showed con- 
siderable differences in thermospore phylotype 
richness, despite an absence of strict endemism of 
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phylotypes that occurred at multiple locations (that 
is, exclusive presence in only one geographic 
region). The number of detected thermospore 
phylotypes per site ranged from 0 to 51 (on average, 
13.3 ± 11.6 (s.d.), 27 = 81) (Supplementary Table Si). 
Thermospore phylotype richness was generally 
higher in Northern Europe (36.4 + 10.2 per site, 
27 = 5), the Guaymas Basin (28.7 + 18.6 per site, 
27 = 3), West Svalbard (17.9 + 8.3 per site, 27 = 23) 
and South Asia (22.8 + 15.2 per site, 27 = 4), whereas 
sites in the Baffin Bay (5.4 + 4.7 per site, 27 = 25), the 
South Pacific (9.5 + 6.5 per site, 27 = 4), East Sval- 
bard (10.3 + 6.2 per site, 27 = 4) and other regions 
(7.2 + 5.7 per site, 27 = 13) showed comparatively 
lower phylotype richness (Figure 1). These differ- 
ences in site occupancy indicate non-random 
variation in the distribution of thermospore 
phylotypes across the investigated sites. Out of 
the 146 thermospore phylotypes, 21 were detected 
at more than 15 locations and were widely 
distributed across different oceanic regions ('cos- 
mopolitan phylotypes', Figure 2, Supplementary 
Figures S4 and 5). This cosmopolitan distribution 



of thermospore phylotypes suggests that these taxa 
(i) are globally dispersed by ocean circulation 
and/or (ii) derive from multiple, globally wide- 
spread source environments that support growth 
of similar communities of thermophilic, endo- 
spore-forming bacteria. The cosmopolitan ther- 
mospore phylotypes were related to Bacillus 
(TSP005, TSP010, TSP013, TSP021), the Aeriba- 
cillus-Geobacillus lineage (TSP003, TSP007, 
TSP016), Desulfonispora (TSP014), Desulfotoma- 
culum (TSP004, TSP006, TSP015), Clostridium 
(TSP018, TSP019, TSP020), the Brassicibacter- 
Sporosalibacterium lineage (TSP002, TSP009, 
TSP012), or were only assigned at the family or 
order level [Clostridiales Family XI. Incertae 
Sedis: TSP008; Christensenellaceae: TSP017; 
Bacillales: TSP001, TSP011) due to lack of close, 
cultivated relatives (Figure 2). In contrast, 82 
thermospore phylotypes were found only at five 
or fewer locations (Supplementary Figures S3 and 
4), indicating more restricted occurrence and/or 
low abundance below the detection limit of our 
approach. 
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Figure 2 Phylogeny and geographic distribution of cosmopolitan thermophilic endospores. 16S rRNA-based phylogenetic tree of 
thermospore phylotypes (TSPs) that were detected at 15 or more sites. Scale bar indicates 1% sequence divergence as inferred from 
RAxML. Phylogenetic affiliations are highlighted in different colors. Colored bars indicate broad geographic regions where the 
thermospore phylotypes were present. Numbers indicate the number of sites at which a thermospore phylotype was detected, whereas 
the height of each bar color indicates the number of those phylotypes detected in each ocean region. TSP, thermospore phylotype. See 
Supplementary Figure S3 for an extended tree showing the phylogeny and geographic distribution of all 146 thermospore phylotypes. 
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Global biogeography patterns of thermophilic 
endospores show the influence of dispersal limitation 
To analyze the relationship between geography and 
differences in endospore community structure, we 
plotted similarity between thermophilic endospore 
communities vs geographic distance for each 
location pair. This analysis revealed a significant, 
negative distance-decay curve (linear regression of 
loglO-transformed variables with slope = —0.68177, 
intercept = -0.09818, P<0.001) (Figure 3). This 
result shows that thermophilic endospore commu- 
nities are non-randomly distributed, which suggests 
that dispersal limitation is influencing the beta- 
diversity of thermospore phylotypes. Alternatively, 
the shape of the distance-decay curve may be 
influenced by differences in the ability of endospore 
phylotypes to resist decay during dispersal (Nicholson 
et al, 2000; McKenney et al, 2013). To account for 
this, we investigated if the distribution of individual 
phylotypes, which should be less biased by within- 
phylotype variances in endospore survival, is 
significantly different from a random geographic 
distribution. Twelve of 146 thermospore phylotypes 
(TSP001, TSP007, TSP008, TSP009, TSP011, TSP015, 
TSP020, TSP027, TSP060, TSP061, TSP096, TSPlll) 
were seemingly dispersal-limited because they had 
significantly geographically clustered occurrence; that 
is, the average geographic distance between the 
locations of each phylotype was significantly less 
than would be expected from its random occurrence at 
the same number of locations. 



Correlation network analyses reveal a widely distributed 
core of frequently co-occurring thermophilic endospores 
Network analysis was recently incorporated 
into microbial biogeography research to explore 
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Figure 3 Distance decay of thermophilic endospore commu- 
nities. The relationship between thermophilic endospore com- 
munity similarity (Jaccard similarity, loglO-transformed) and 
geographic distance (in kilometers, loglO-transformed) is shown. 
The linear regression of log-transformed variables is significant 
[P< 0.001) and has a negative slope [m = — 0.68) and shows that 
communities of thermophilic endospores are non-randomly 
distributed. 
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co-occurrence patterns of microbial taxa (Fuhrman, 
2009; Barberan et al, 2012). We used correlation 
network analysis to better describe the biogeography 
patterns of thermophilic endospores identified by 
the distance-decay curve and investigated which 
phylotypes tend to co-occur; that is, which phylo- 
types are found together at some sites and are 
commonly absent from others. The co-occurrence of 
thermospore phylotypes was calculated using non- 
parametric Spearman correlations of thermospore 
phylotype presence/absence across all sampling 
sites and, to explore groups of co-occurring phylo- 
types, network analysis was used for visualization. 
We identified eight small networks of different 
complexity, that is, numbers of nodes and edges 
(connections) (Supplementary Figure S6A). The 
largest network was characterized by four highly 
connected, central phylotypes (^4 connections; 
TSP004, TSP008, TSP009, TSP015) and nine per- 
ipheral phylotypes that occurred at several locations 
(Figure 4a). Once we identified the most frequently 
co-occurring phylotypes, we explored similarities 
between thermospore phylotype communities by 
inferring association networks of sediment loca- 
tions. The largest network contained 18 locations, of 
which many harbored several thermospore phylo- 
types (Figure 4b), whereas 10 other, much smaller 
networks were composed of ^4 locations 
(Supplementary Figure S6B). The 18-location net- 
work was composed of most West Svalbard samples, 
but also samples that were geographically very 
distant from the Svalbard archipelago, namely two 
samples from the Baltic Sea (Aarhus Bay, Arkona 
Basin) and one sample from the Gulf of California 
(Guaymas Basin) (Figure 4b). This network contains 
a phylogenetically diverse core of co-occurring 
thermophilic endospores that are widely but not 
ubiquitously distributed in the oceans (for example, 
one of the phylotypes, TSP001, was detected in 39 
sediments). This raises questions about whether 
certain thermospore phylotypes share common 
source environments from where they are dissemi- 
nated in a non-random manner over long distances 
and along similar, as yet unidentified travel routes. 



Hydrothermal sediments of the Guaymas Basin as 
potential source environments for endospore- forming 
thermophiles 

Based on the identity, physiology and sedimentation 
rates of thermophilic spores in Svalbard and Aarhus 
Bay, we have previously postulated biogeochemical 
and geological characteristics of hypothetical source 
environments for these bacteria (Hubert et ah, 2009, 
2010; de Rezende et ah, 2013). First, such anoxic 
marine environments must be warm enough to allow 
vegetative growth of the diverse community of 
anaerobic, endospore-forming thermophiles. Second, 
sufficiently strong fluxes (e.g., fluid transport) from 
these environments into the water column must 
physically transport cells into circulating seawater. 
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Figure 4 Network analysis of thermophilic endospore co-occurrence (a) and location (b). (a) Largest network of co-occurring thermospore 
phylotypes. Each node represents a thermospore phylotype (TSP). Presence of an edge between two nodes shows a strong correlation 
between these two phylotypes, which is indicative for co-occurrence. Circle size indicates site occupancy Map shows locations where 
phylotypes of the network are present. The circle color indicates how many of the 13 phylotypes comprising the network are present at this 
site, (b) Largest location network. Each node represents a location. Presence of an edge between two nodes corresponds to a high Bray-Curtis 
similarity (^0.6) between the endospore communities at these two locations. Circle size indicates thermospore phylotype richness. Maps 
show the global locations and correlations of these sites. The map on the right shows a magnification of Svalbard for enhanced resolution. 



Consistent with this general description are, for 
example, pressurized gas or oil reservoirs in seabed 
sediments, and parts of oceanic spreading centers 
including hydrothermal vents and the sediments 
overlying them (Head et al, 2003; Cowen, 2004; 
Brazelton et al, 2006; Judd and Hovland, 2007; 
Hutnak et al, 2008; Hubert and Judd, 2010). 
Anaerobic thermophiles in these hot sediment and 
ocean crust habitats are connected with the 
overlying cold ocean water via advection of seabed 
fluids such as hydrocarbons and recirculating water 
discharged from the ocean crust (Judd and Hovland, 
2007; Hutnak et al, 2008). 

The hydrothermal surface sediments and the 
general geology and geography of Guaymas Basin 
satisfy these criteria (Supplementary Results and 
Discussion). Accordingly, a very large proportion 



(44.5%) of all 146 thermospore phylotypes were 
detected overall in the three Guaymas Basin 
sediments investigated. Although surveys with 
much finer genetic resolution than offered by the 
16S rRNA gene approach are required to prove that 
identical microorganisms co-occur at distant loca- 
tions, it is intriguing that the Guaymas Basin 
sediment community comprised 15 of 25 cosmopo- 
litan thermospore phylotypes and 8 of the 13 
co-occurring phylotypes (including all four highly 
connected phylotypes in the largest network) 
(Figure 4a). High thermospore phylotype richness 
despite restricted access due to the geographic 
isolation of the semi-enclosed Gulf of California 
supports the hypothesis that environments in the 
Guaymas Basin are sources of thermophilic spores 
to the overlying water column. Thermohaline water 
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circulation in the Gulf of California is primarily 
influenced by wind, which alters directions and 
water transport between the gulf and the Pacific 
Ocean seasonally (Bray, 1988). It is conceivable that 
Guaymas Basin-derived spores embark on and 
survive long-distance travel to far away Svalbard 
and North Sea sediments via the global thermoha- 
line circulation (Rahmstorf, 2000) (Figure 1). How- 
ever, it may be more likely that these cold 
destinations are supplied with thermophilic spores 
from nearby warm environments that were not 
investigated in this study, and that dispersal barriers 
may prevent thermophilic endospores that escape 
the hot anoxic sediments in the Guaymas Basin from 
reaching far away water masses in the north 
Atlantic. Possible examples of potential source 
environments closer to Svalbard and the Baltic Sea 
could include systems such as the Lost City 
hydrothermal vent field associated with the mid- 
Atlantic ridge (Kelley et al., 2005; Brazelton et al., 
2006) or the oil-bearing sediments in the North Sea 
(Gittel et al., 2009), which both support thermo- 
philic sulfate-reducing Desulfotomaculum species 
that are closely related to the thermospore phylo- 
types identified in this study. It is conceivable that 
thermophilic bacteria get distributed along the mid- 
oceanic ridges throughout the world oceans, as 
exemplified by hydrothermal vent animals and their 
symbionts (Petersen et al, 2010; Vrijenhoek, 2010; 
Teixeira et al., 2011), creating an extended series of 
source habitats from where thermophilic endo- 
spores may be further dispersed to the cold 
sediments tested in this study. The location network 
shown in Figure 4b, thus likely represents a fraction 
of a much larger network of warm sources and warm 
or cold destinations for dispersing thermophilic 
endospores that are interconnected through ocean 
currents. 

Anthropogenic sources may also contribute to 
the presence of anaerobic thermophiles in marine 
sediments. North Sea oil production systems have 
been proposed as possible sources for thermophilic 
Desulfotomaculum species in Aarhus Bay sediments 
(de Rezende et al., 2013) and close relatives to 



thermospore phylotypes detected in the present 
study were found in oil production water (Gittel 
et al., 2009). Close relatives were also detected in 
several anthropogenic, as well as natural terrestrial 
habitats such as bioreactors, activated sludge, com- 
post and soil environments (Supplementary Tables 
S3 and 4). However, the phylotypes we have 
identified were enriched in artificial seawater, 
which is suggestive of a marine origin. Also, the 
presence of a widely distributed thermophilic Desul- 
fotomaculum alkaliphilum phylotype (Figure 2) at 
depths corresponding to ~4500 years of sedimen- 
tation in Aarhus Bay (de Rezende et al., 2013) 
indicates that there is a natural long-term dispersal 
of thermophilic endospores occurring in the marine 
environment that is independent of human activity. 



Differences in thermophilic spore richness and 
phylogeny in Arctic regions correlate with connectivity 
of major water masses 

The variable compositions of thermophilic endo- 
spore communities in Arctic sediments from the 
Baffin Bay and the east and west sides of the 
Svalbard archipelago suggest that thermospore phy- 
lotypes are not equally dispersed across these Arctic 
regions. Sediments from the Baffin Bay contained a 
significantly lower number of thermospore phylo- 
types than sediments from West Svalbard (Table 1). 
While sediments from only four locations on the 
East coast of Svalbard were available, we detected a 
similar yet not statistically significant trend of lower 
richness in East compared with West Svalbard 
sediments. Analysis of more samples from East 
Svalbard is required to confirm this trend. We 
additionally revealed a contrasting geographical 
distribution in the potential for thermophilic 
sulfate reduction that was generally consistent with 
the observed differences in thermospore phylotype 
richness (Figure 1). Thermophilic sulfate reduction 
was significantly more prevalent in West Svalbard than 
in East Svalbard and Baffin Bay sediments (Table 1). 

The differences in thermophilic spore richness 
and metabolic potential were also reflected by the 



Table 1 Comparison of thermospore phylotype richness, phylogenetic composition and prevalence of thermophilic sulfate reduction 
between Arctic regions 



West Svalbard East Svalbard Baffin Bay West Svalbard vs West Svalbard vs East Svalbard vs 
(n = 23) (n = 4) (n = 25) Baffin Bay** East Svalbard 11 Baffin Bay" 



Thermospore phylotype richness 0 17.9 ±8.3 10.3 ±6.2 5.4 ±4.7 P<0.001 Not significant Not significant 

Clostridiaceae richness 6.0 ±2.8 2.0 ±1.8 1.5 ±2.1 P<0.001 P=0.02 Not significant 

Peptococcaceae richness 3.4 ±2.9 2.0 ±1.6 0.6 ±0.9 P<0.001 Not significant Not significant 

Bacillaceae richness 2.2 ±1.9 3.0 ±1.4 2.4 ±1.4 Not significant Not significant Not significant 

Thermophilic sulfate reduction 87 25 24 P<0.001 P=0.02 Not significant 

(% of sites tested positive) 



""Significant differences between Arctic regions were determined by pairwise comparisons using the Mann-Whitney U-test for richness and 
Fisher's exact test for thermophilic sulfate reduction. P-values ^0.05 were considered significant. 

b All significant P-values in this column are also significant (P<0.01) when only surface (0-10 cm, u = 13) or only deep (10-121 cm, u = 12) 

sediment samples are analyzed. 

C A11 richness values are presented as mean±s.d. 
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phylogenetic composition of thermophilic endo- 
spore communities in these Arctic regions (notwith- 
standing the aforementioned weak indications for 
strict endemism for any individual phylotypes). 
Of the three major families, West Svalbard sediments 
exhibit significantly more thermospore phylotypes 
belonging to the Clostridiaceae and the Peptococcaceae 
than Baffin Bay sediments, whereas Bacillaceae 
phylotypes were not differentially distributed between 
the Arctic regions (Table 1). Again, we observed a 
similar pattern between West Svalbard and East 
Svalbard, although differences were only significant 
for Clostridiaceae phylotypes (Table 1). No significant 
differences between East Svalbard and Baffin Bay 
sediments were found among these three taxonomic 
families (Table 1). Furthermore, frequently co-occurring 
thermospore phylotypes (e.g. TSP004, TSP008, 
TSP009, TSP015), as identified by network analysis 
(Figure 4), were almost exclusively limited to West 
Svalbard sediments. The observed distribution pattern 
is not impacted by the fact that 12 of 25 Baffin Bay 
sediment samples originate from a sediment depth 
below 10 cm (Supplementary Table Si). Restricting our 
analyses to only the shallow and deeper Baffin Bay 
sediments showed (i) no significant difference between 
them (e.g., thermospore phylotype richness: surface 
samples = 4.8 ±4.7, deeper samples = 6.1 ±4.9) and 
(ii) that both had the same significant differences 
compared with West Svalbard regarding richness, 
phylogenetic composition and prevalence of thermo- 
philic sulfate reduction (Table 1). 

One dispersal-limiting factor that might be 
responsible for the different thermophilic endospore 
communities in Svalbard and Baffin Bay sediments 
is the rate at which spores are deposited to the 
seafloor. Based on sediment thickness and age, 
sedimentation rates in the Baffin Bay and off the 
coast of Svalbard are estimated to be in the same 
order of magnitude (5-25 cmka" 1 , Kallmeyer et ah 
(2012)). In contrast, reports for a few locations in 
these regions suggest that sedimentation rates are 
lower in the Baffin Bay (6.5 cmka -1 , Simon et al. 
(2012)) than in Svalbard fjords (180 cmka" 1 , Hald 
et al (2001); 190 cmka" 1 , Hubert et al (2009)). 
However, it is currently unknown how efficient 
endospores in the water column are deposited to the 
sediment, that is, how the specific rate of spore 
sedimentation relates to the overall sedimentation 
rate. It is noteworthy that for most locations our 
enrichment inoculum derives from homogenized 
sediment from 0 to 10 cm depth and therefore 
integrates endospore communities that accumulated 
over 50 to 1500 years of deposition. Given the 
longevity of thermophilic spores (de Rezende et al., 
2013), our approach may diminish the impact of 
sedimentation rates on the recovered richness of 
thermophilic spore communities. 

Besides the impact of local hydrography, sedi- 
mentation, and the possibility that unknown, local 
sources contribute to biogeography of thermophilic 
endospores, we hypothesized that dispersal-driven 



distribution also depends on how local currents in 
the sampled regions are connected to global ocean 
circulation. The observed West-East difference in 
thermophilic endospore community structure in 
Svalbard appears to correlate with the regional 
hydrography of the archipelago, although the differ- 
ence is not statistically significant in all parameters 
(Table 1). Two major ocean currents define the 
hydrography of Svalbard. The warm West Spitsber- 
gen Current, which is an extension of the Gulf 
Stream, flows from the Atlantic Ocean northwards, 
whereas the East Spitsbergen Current transports 
cold water from the Arctic Ocean southwards 
(Loeng, 1991) (Figure 1). Although the core of the 
West Spitsbergen Current that flows along the 
continental slope is separated from the West 
Spitsbergen shelf waters by the Arctic Front, 
extensive cross-front exchange takes place below 
50 m depth (Saloranta and Svendsen, 2001) and it 
has been shown that the West Spitsbergen Current 
extends its influence deep into the fjord system of 
West Svalbard (MacLachlan et al, 2007). 

The formation of deep and bottom water masses in 
the Baffin Bay are not fully understood (Tang et al., 
2004), but they are generally influenced by cold 
Arctic Ocean water via the Baffin Island Current in 
the north and by the Atlantic Ocean via the West 
Greenland Current (Figure 1). The low thermospore 
phylotype richness in the Arctic water-influenced 
sediments in the Baffin Bay suggests a lower influx 
of thermophiles derived from the Arctic Ocean. 
Reduced thermospore phylotype richness in the 
Arctic Ocean is possibly due to limited thermophile 
sources and lower connectivity to water masses from 
the rest of the world oceans where thermophile 
sources may be far more abundant. 

Recent studies of vegetative microbial commu- 
nities in the global ocean have suggested that 
hydrography and geographic isolation of the Arctic 
Ocean represents an effective dispersal barrier for 
microorganisms (Galand et al., 2010; Ghiglione 
et al., 2012; Hamdan et al, 2013; Sul et al, 2013). 
By tracking thermophilic endospores whose Arctic 
biogeography is only controlled by passive disper- 
sal, our results demonstrate this dispersal limitation 
to be true. Importantly, evidence presented here 
suggests that this dispersal limitation even holds 
true for marine bacteria with enhanced survival 
capacities that are less prone to death during long- 
term and long-distance dispersal under unfavorable 
environmental conditions, that is, perhaps the most 
likely candidates to not be dispersal-limited. Different 
phylotypes of globally widespread thermophilic 
endospores co-occur at multiple distantly separated 
locations and thus could be simultaneously 
dispersed from common sources and/or along 
similar paths according to non-random mechanisms. 
Global dispersal routes and frequencies of marine 
microbial taxa are largely dependent on the degree of 
connectivity between major oceanic water masses. 
Dispersal barriers of this kind exert greater control on 
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the composition of marine microbial communities in 
water masses that participate minimally in global 
ocean circulation such as the Arctic Ocean. 
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